Quantum Cloning for Absolute Radiometry 
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In the quantum regime information can be copied with only a finite fidelity. This fidelity gradually 
increases to 1 as the system becomes classical. In this article we show how this fact can be used to directly 
measure the amount of radiated power. We demonstrate how these principles could be used to build a 
practical primary standard. 



o 



a\ 



I 



> 

in 

00 

en 

o 
o 

> 

X 



Since its inception quantum mechanics has had a deep 
tie with radiometry, the science of measurement of electro- 
magnetic radiation. The electrical substitution radiometer, 
developed by Lummer and Kurlbaum in 1892 fT], was used 
to observe the spectral distribution of a heated black body. 
In 1900 Max Planck was able to describe this distribution 
by assuming that electromagnetic radiation could only be 
emitted in multiples of an energy quantum E = hu. This 
discovery not only provided an accurate law relating the 
radiated spectral density to temperature, but laid the foun- 
dations of quantum physics. The electrical substitution ra- 
diometer is still used as the primary standard for spectral 
radiance by many metrology laboratories. These systems 
have been improved over more than a century and can now 
achieve absolute uncertainties better than 10""', when op- 
erated at relatively high powers [3- 

More recently, nonlinear optical effects such as Sponta- 
neous Parametric Down Conversion have provided a new 
primary standard based on the correlations of quantum 
fields fSj. The accuracy of these techniques has improved 
by nearly one order of magnitude every ten years, and is 
currently of the order of 10^'^. These systems are limited to 
the photon-counting regime, with a recent theoretical pro- 
posal for extension to higher photon rates [4J. 

In this letter, we present a radiometer that overcomes 
these limitations and works over a broad range of pow- 
ers: from the single photon level, up to several tens of nW 
(~ 10^^ photons/s), i.e. from the quantum regime to the 
classical regime. In fact, our system is able to provide an 
absolute measure of spectral radiance by relying on a par- 
ticular aspect of the quantum to classical transition: as the 
number of information carriers (photons) grows, so does 
the fidelity with which they can be cloned. For an opti- 
mal cloning machine Ii5rt9j this relation can be derived ab 
initio lITOl fTTl so that a measurement of the fidelity of the 
cloning process is equivalent, as we shall see below, to an 
absolute measurement of spectral radiance. 

Optimal cloning has been demonstrated in a variety of 
systems ll6l[8j|9l. Stimulated emission in atomic systems is 
particularly practical as high gains can be easily achieved 
and the entire system can be implemented in-fibre which 
both ensures the presence of a single spatial mode and 
makes the system readily applicable, though not limited, 
to current telecom technology. 

Principle of operation - The aim of this experiment is to 
produce an absolute measurement of luminous power Pin. 



We will do this by using an optimal Universal Quantum 
Cloning Machine (QCM). As we shall see such a device is 
able to directly relate a relative measurement of two or- 
thogonal polarizations at its output to P{^. The relative 
measure that we use is the fidelity T which is the mean 
overlap between the input and output polarization, and can 
be expressed as follows: 

where P\\ and P±^ are the output powers in the polarizations 
parallel and perpendicular to the polarization of the input 
light. 

For an optimal QCM (losses are considered in the next 
section), the fidelity of a cloning process from A^ to M 
qubits can be derived ab initio 1,10,1 to be: 

_ NM + N + M 
-^^^*' " NM + 2M • ^^^ 

This equation remains valid when we clone a large num- 
ber of polarization qubits distributed over a large number 
of temporal modes and can be rewritten in terms of the av- 
erage number of input and output photons per (temporal) 
mode^Uin and /Xout 111: 

_ A^ll _ /^in/J-out + /^in + /^-out 



'J II. 



(3) 



/-t|| + /^_L Min/^out + 2^out 

where ^Uout contains both a number of exact copies of the 
input signal and intrinsic noise due to the amplification pro- 
cess, i.e. /iout = ^11 + Ai±- 

It is also possible to express /Xout as a function of /Xin and 
the amplifier gain G [ 12 1 : /^out is the sum of the stimulated 
emission G fi^ and the spontaneous emission, equivalent 
to amplifying the vacuum, so that: 

Mout = G/Xi„ + 2(G-l). (4) 

Equations Q and Q can be combined to obtain the spec- 
tral radiance //;„ as a function of fidelity and gain: 

_2J='G-G-2J='+1 2F-1 

^"^ ~ G-FG - Y^y ^^^ 

with the approximation holding for G ^ 1. P\a can be de- 
rived from /Xin and a measurement of the number of modes 
per unit time r~^. 

Three aspects make this scheme attractive: the first is 
that after amplification input power information is polar- 
ization encoded and is therefore insensitive to losses [131; 
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FIG. 1. Experimental arrangement: a broadband source with controllable polarization and power is amplified by an Erbium doped 
fiber amplifier. The degree of polarization (DOP) of the amplified light is then measured with a polarimeter. Spectral bandwidth is 
determined by a tunable filter. A value for the input power can be calculated from the DOP and compared with a calibrated powermeter 
(monitor). 



the second is that the experiment can be performed entirely 
in fiber, ensuring the selection of a single spatial mode. The 
third advantage is that this scheme works over a broad scale 
of powers: from single photon levels up to several tens of 
nW (~ 10" photons/s). 

Non-ideal cloning - The reasoning presented above as- 
sumes the universal cloning process to be optimal. It has 
been shown theoretically that amplification in an inverted 
atomic medium indeed provides optimal cloning [5|, but 
for precision applications it is important to consider the 
possible effects of a non perfectly inverted medium, which 
we model by a succession of infinitesimal gain and loss 
elements, G„ and ?7„, as shown in Fig. [Sj^a). We have 
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FIG. 2. (a) Model of a non-totally inverted medium as succession 
of infinitesimal gain elements Gn and loss elements r\n. (b) Each 
loss element r\n within the fiber is equivalent to a smaller loss 
element Qn before the amplifier. 



shown lfT4ll that this model is equivalent to a much sim- 
pler one. Each loss element ?]„ can be represented by a 
different loss element Q^ before an optimal cloning ma- 
chine with gain G' , as shown in Fig.jlj^b). It can be shown 
that the product Q of all Q„ can be expressed as: 



^/ II ^n ) ^n 



GoVn 



G^r]^ + (1 - r?„ 



(6) 



where Gq is the effective gain between the beginning of the 
amplifier and element ??„. A fully inverted medium would 
have Q = 1. From Eqn. (|6]l it is apparent that the effect of a 
small loss {rfn ^ 1) is proportional to l/Gp . As G^ grows 
exponentially over the length of the fibre, losses towards 
the end can be neglected. At the beginning of the ampli- 
fier two effects guarantee that the medium is fully inverted: 
the input signal is small, as it has not been amplified yet, 
and the signal and pump co-propagate, ensuring maximum 
pump power in this region. Cloning optimality can then be 
achieved in a non-ideal amplifier. 

Experimental arrangement - Fig. [T] shows the setup, 
which can be conceptually divided in three main parts: 
generation of a set amount of power, amplification and fi- 
delity measurement. To test our system, we prepare states 
with a known number of photons per mode (^in). This is 
done using a polarized LED that is passed through a po- 
larization controller (scrambler) and a variable attenuator. 
The power is then split (50:50), with one branch moni- 
tored on a calibrated powermeter, while the other is sent to 
the amplification stage. Amplification is provided by 2 ni 
of Er^ doped fiber, (attenuation 16.7dB/ni at 1530 nin), 
pumped by a 980 nm diode laser. The pump light is com- 
bined with the signal on the input of the Er'^^fibre using 
a wavelength division multiplexer (WDM), and an isolator 
is placed before the input to prevent unwanted resonances. 
After the Er'^^doped fibre, most of the pump power is re- 
moved using an additional WDM. In this realization, the 
no-cloning theorem is guaranteed by the Er^^ spontaneous 
emission, which adds randomly polarized photons to the 
signal. We used an optical frequency-domain reflectome- 
ter lITSl to verify that the gain per unit length is constant 
over the entire fibre, indicating that the atomic medium 
is fully inverted. Results are shown in Fig. [3] The mea- 
surement stage consists of a grating-based tunable filter 
and a polarimeter. The filter has a width of 273.3(5) pm 
(FWHM), which ensures that the effects of polarization 
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FIG. 3. Optical frequency-domain reflectometer measurement 
showing homogenous gain per unit length within the Er'^^ doped 
fibre. The solid line is an exponential fit of the data. 



mode dispersion can be neglected. The polarimeter mea- 
sures the degree of polarization (DOP) with a nominal ac- 
curacy of 1%, where the DOP is defined as the polarized 
power (in any basis) Ppoi over the total power Ptot, and is 
related to fidelity by J" = (1 + D0P)/2. Using a po- 
larimeter rather than simply a polarizing beamsplitter and 
powermeter is less accurate, but allows us to test whether 
the system works equally well for arbitrary input states of 
polarization, i.e. whether the QCM is truly universal. 

Experimental procedure — To evaluate the accuracy of 
our system, we will need to compare our measurement of 
/Xin with the value ^l^ obtained from the reference power- 
meter. To do so, we first measure the ratio between the 
power at the monitor output and the power at the entrance 
of the amplifier within the bandwidth of the tunable filter. 
This is done by placing the filter just before the amplifica- 
tion stage (position 'a' in Fig.[T]l. Together with a measure- 
ment of the filter's attenuation and bandwidth, this allows 
us to obtain ^l^ from the monitor power. The filter is then 
placed after the amplification stage (position 'b') so that 
Er'^^ spontaneous emission outside of the bandwidth of in- 
terest is eliminated. We then vary ^l^ using the attenuator, 
and record the monitor power versus the degree of polar- 
ization. For each jjl*^ the measurements are repeated for 20 
different input polarizations, to estimate uncertainties. 

Results -The fidelity of the cloning process is a measure 
of spectral radiance. In order to measure power it is neces- 
sary to have an accurate measure of the number of modes 
involved. Using a single-mode fibre ensures that there is 
only a single spatial mode: only the number of temporal 
modes per second need to be measured. It is convenient to 
define the coherence time r^ as in ||T6l : 



\l{r)\'dT 



(7) 



where 7(r) is the autocorrelation function normalized such 
that 7(0) = 1. Using this definition, the coherence length 



c Tc is the length of the unit cell of photon phase space [T^, 
so that the number of modes per second is simply t^^- 
Measuring this value with an optical low-coherence inter- 
ferometer (Fig. [4I1 yields r^ = 19.71(4) ps which corre- 
sponds, assuming a Gaussian shape, to wavelength FWHM 
of AA = 273.3(5) pm. We also performed a (less precise) 
spectrometric measurement yielding AA = 271 pm. With 
this filter, a mean of one photon per temporal mode corre- 
sponds to 6.461 nW. 




FIG. 4. Autocorrelation function of the source after the filter. 
7(r) is the fringe visibility measured with a low-coherence inter- 
ferometer; Tc will simply be the numerical integral of this data. 
The solid line is a Gaussian fit. 




FIG. 5. Fidelity versus number of input photons per mode, fitted 
with Eqn. iBl. Representative errors are shown as boxes on some 
points. The inset shows the output versus input number of pho- 
tons per mode, the line is a fit on the first data points (/xj^ < 1), 
showing reduced gain as /ijn grows. 



We measure the amplifier gain by directly comparing the 
power at the output of the amplifier with the power at the 
input. The inset of Fig. [5] shows a typical plot, in terms 
of iil^ and /i*m, the thickness of the line represents random 
errors. Note that the gain is G = dfiout/dfim', so that any 



systematic error in either the power measurement or the es- 
timation of the number of modes cancels. The line in the 
inset of Fig. Blis a fit of the data for /ij*„ < 1, revealing 
that at high fii^ the gain is reduced. This effect could be 
minimized by pumping from both sides of the Er'^^doped 
fiber. Nevertheless, the gain is constant for /Uin < 2, al- 
lowing us to assume within this range that the intercept fiQ 
corresponds to the spontaneous emission (2 G — 2) from 
Eqn. (J4|, so that: /Xout = G fXi„ + fiQ. In this range it is 
then possible to measure fj,[„ without distinguishing the po- 
larizations, as Hin = 2(/i*m - fJ-*n)/l^o - 2. 

We then measure the fidelity J^ versus /J,*^; Fig. |5] shows 
a typical plot, which can be fitted with Eqn. Q, where ^in 
has been replaced with kfi*^, and k is the fitted parameter. 
With this definition k = //in//i*, represents the discrepancy 
between our measurement of /Xin and the value fj,'^^ obtained 
from the reference powermeter. Here, k also accounts for 
the possibility of non-optimal cloning which would intro- 
duce a further factor Q < 1 equivalent to a loss on the 
input of the cloning machine. The fitted curve in Fig. [5] 
yields k = 1.013(5), where the error indicated represents 
statistical uncertainty. Systematic errors, as we shall see 
in the next section, could be up to one order of magnitude 
higher. 

Error estimation - The aim of this experiment was to 
demonstrate the principle of a cloning radiometer, rather 
that to build a standard that can compete with metrology 
laboratories. It is however important to discuss the errors 
involved, both for the interpretation of the results and to 
evaluate the applicability of this method. One of the ad- 
vantages of this technique is that relative measurements, 
which usually have small errors, are used; but how does a 
small uncertainty in the fidelity A J-" translate into an error 
in the measurement of ^in? From Eqn. Q we obtain: 



A/iin = (2 + /i,n)'A^, 



(8) 



so that A^in//iin has a minimum of Afiin/fim = 8AJ^ at 
/Xin = 2, i.e. when spontaneous and stimulated emissions 
are equal. At higher spectral radiances, Ai^iin/nm rises lin- 
early with //in- The spectral bandwidth of the filter can 
be chosen as to operate in the desired power regime: our 
system is optimal at 13 nW, commercially available filters 
would allow this point to be easily lowered to 100 pW. 
From preliminary tests we estimate that this technique 
would work to an upper limit of 100 nW, above which the 
effects of polarization mode dispersion and wavelength de- 
pendence of the components need to be taken into account. 
The two main systematic uncertainties in our system are 
due to the reference powermeter, and to the polarization 
measurement. The powermeter is an EXFO PM-1100, re- 
cently calibrated by METAS to an absolute uncertainty of 
0.7% and with a measurement to measurement standard de- 
viation of 0.5% (including fibre re-connection). The linear- 
ity of this powermeter is within this uncertainty over its en- 
tire range. The fidelity is measured with a Profile PAT 9000 



polarimeter which has a nominal AJ^ = 0.5%. We no- 
ticed that the fidelity was overestimated by 1 % for unpolar- 
ized light, and underestimated by 0.2% for polarized light. 
Systematic errors are dominated by the polarimeter, so that 

A/iin//Xin ~ 4% for flin = 2. 

Conclusion - We have shown that the fidelity of 
cloning can be used to produce an absolute power mea- 
surement with an uncertainty only limited by the uncer- 
tainty of a relative power measurement. We demonstrate 
the scheme with an all-fiber experiment at telecom wave- 
lengths, achieving an accuracy of 4%, with much space for 
improvement by a metrology laboratory. The experiment 
also demonstrates optimal 1 — )• 67 cloning, and is an in- 
teresting application of Quantum Information Science, and 
in particular of the study of the quantum to classical transi- 
tion. 
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Appendix A: treatment of potential losses in a doped fiber 
amplifier 

In order to evaluate the feasibility of optimal universal 
quantum cloning via stimulated emission in an Er'^^ doped 
fiber, we should take into account the potential effects of 
internal losses. The amplification medium can be naively 
modeled, as shown in Fig. [6] a sequence of thin amplify- 
ing atomic slices spaced out by beam splitters, representing 
the internal optical losses. The propagation of the photonic 
mode in a lossy amplifier can be seen as successive inter- 
action with these elements. 
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FIG. 6. Symbolical representation of the lossy amplification of 
mode a as a sequence of distinct amplifying elements (rectangles) 
spaced out by beam splitters (/ signs). \Bn) and \Cn) stand for 
the initial states of the local auxiliary modes. 



The interaction of the input propagation mode a with the 
71-th beam-splitter can be represented by the hamiltonian 
H^ = iX(acl^ — a^Cn), where 'L' stands for losses, A is a 
constant and c„ is an auxiliary mode initially in the vacuum 
state (|C„) = |0)^ ). By using the time evolution operator 



c/i 



in the Heisenberg picture the action of the 



n-th beam-splitter on modes a and c„ can be expressed by 
the relation 



u!:^au^ 
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+ yr 



'jn C-n 5 



(9) 



where ?]„ is the specific transmission coefficient of the 
beam-splitter element. 

A similar relation can be found for the amplifying ele- 
ment. It has been shown in ifTTl that amplification in an in- 
verted atomic medium provides optimal universal cloning, 
equivalent to stimulated parametric down conversion in 
nonlinear crystals. For this reason the interaction of the 
propagation mode a with an amplifying atomic element 
can be expressed more conveniently by the hamiltonian of 
the parametric amplification process, which for the n-th 
amplifying element is H^ = ixio- bn — a^bl^), where x is a 
constant and 'A' stands for 'amplification'. 6„ in the para- 
metric case is the mode of the anticlones [5|, while in the 
atomic case it represents a collective "desexcitation" of the 



atoms in the nth-slice [17|. Initially \B„ 



|0) 



b,^ 



mean- 



ing that the population is inverted in the element (all the 
atoms are in the excited state). The amplifying elements 
are characterized by different gain values G„. Notice that 
the hamiltonian H;^ has a second term containing aj_ that 
ensures universality, however this term is decoupled and 
doesn't affect mode a. 



By using the time evolution operator U^ 



e 1 , m 

the Heisenberg picture the action of the amplifying element 
on modes a and 6„ yields: 



a^ = U^^aU^ 



(10) 



Now let us consider the two different situations represented 
schematically in Fig. [7j In the first case the propagation 
mode ttin is through an amplifying element of gain G be- 
fore interacting with mode c in a beam splitter of transmis- 
sion r]. In the second situation the order is inverted with the 
parameters G' and 77'. Using (J9]) and ( lOi the value of a at 
the output for the two cases is: 



^out — \/Gr\ ain + 



rJiG-l)b^ + ^/T^c, (11) 



VG' 



'??' a,. 



+ \{G' 



1)6^ + a/G'(1-7?')c. 
(12) 



Suppose that we fix the value of G and r/ and solve for the 
value of G' and r/' that would give the same output a\^^ = 
al^f. The following three conditions must be satisfied 

= GY, 

= (G"-l), (13) 

= G'(l-r/'). 

Since the first equation in ([131 is the difference of the 
other two, the system always has the solution 




G' 

V' 



r/(G-l) + l, 



(14) 



Gr,+ (1-J))- 

It is easy to verify that this solution satisfies G' > 1 and 
< 77' < 1 for any given rj and G, while the inverse is 
not true. The last equation in ( [T3| ) clearly stipulates that the 
condition G'{l—r]') < 1 must be satisfied if we want to re- 
arrange the elements (else it would imply a negative trans- 
mission 77). The consequence of this result is that we can 
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FIG. 7. The situations where the loss comes before or after the 
amplification are physically equivalent if the scalar parameters 

(G, 77) and (G', 77') satisfy ^ 

pull all the beam-splitter elements in Fig. 1 on the left if we 
take care of correctly modifying the characteristic parame- 
ter for each element. It is well known that a combination of 
beam-splitter interactions with modes q is equivalent to an 
interaction with a single mode c being a linear combination 
of Ci with the resulting transmission rate fj = JJ- Vi' the 
same is valid for a series of amplification layers implying b 



and G = YlGi. So the initial process represented in Fig. 1 
can be equivalently seen as a transmission loss fj' before an 
amplification G', or the other way around with G and fj if 
G'{1 — fj') < 1 is satisfied. Another point is that the lossy 
elements at the beginning of the propagation line will con- 
tribute more to the resulting effective loss that those at the 



end, naively one can see from ( 14 1 that pulling a small loss 
element r/„ = 1 — e on the left through Gi, G2, • • • , G„ 
results in an effective loss of: 

Gi . . . GnTIn 



Vn 



Gi...G„?7„ + (1 -7]„ 



1 



Gi . . . G„ 



(15) 



which illustrates the fact that in real experiments it is 
much better to have a copropagating pump laser. To 
be more quantitative we can introduce continuous fields 
x{z) and X{z), in such a way that propagating from z 
to z + 5z the mode a undergoes an infinitesimal gain 
5G{z) = 1 + xiz)6z and an infinitesimal transmission 
loss 6ri{z) = 1 — X{z)6z (for two infinitesimal elements 
the order is not important) as illustrated in the first line of 
the Fig. [8] Then starting from the left we put all the loss 
elements before the gain elements. We then get a resulting 
equivalent loss at point z of rj{z) and a resulting equivalent 
gain G{z), that would be the second line of the Fig.^ Then 
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FIG. 8. The infinitesimal loss and gain elements are successively 
absorbed in a global gain parameter G{z) and a global loss pa- 
rameter r]{z). In the limit Sz ^ Eqns. ( |14| i yield differential 
equations for these functions. 



using the equations (14i we permutate G{z) with 5ri{z) 
obtaining G'{z) and oij{z); and finally combine them into 



r/(z + 5z) = r]{z)5r]'(z) and G{z + 5z) 
(last line in Fig. IHll to obtain: 



G(z + 6z) 
T]{z + 5z) 
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{dr]{z){G{z)-l) + l)5G{z) ' ' 



In the limit 5z — )• expanding the expressions in the r.h.s. 
to the first order in 6z we obtain the system of differential 
equation that G{z) and rj{z) obey. 



5.G(z) 
d,r]{z) 



ixiz) - X{z))G{z) + Xiz) 
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(17) 



G{z) 



the solution of this system is: 



rj{z) = (1 + /; A(z°)e-i'o^'(x(^^)-^(^^»'^^^)-i 

This solution can be used to take into account the effect of 
any arbitrary loss profile. For our current system the ef- 
fect of losses is negligible, however it is important to know 
that for future systems the loss characteristics of the fiber 
and amplification process can be fully taken into account 
by measuring only two parameters of the fiber: G and r], 
and will thus not impose a limit on the accuracy of a high 
precision system. 

To illustrate the solution of this equation we consider a 
sample with three different profiles of atomic inversion. 



namely Xi{z) 



„~\-a.z 



X2 = cte and Xsiz) 



2 — e"^^ ^'. The total equivalent loss rj for each of these 
three profiles is given in Fig. 4. 




max (loss /gain) 

FIG. 9. Global resulting loss rj for, from top to bottom, xi, X2 
and X3 as a function of a. The total equivalent gain for these 
plots is approximately 50. 



